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Abstract 

The  performance  of  solid  oxide  fuel  cells  (SOFCs)  is  affected  by  various  polarization  losses,  namely,  ohmic  polarization,  activation 
polarization  and  concentration  polarization.  Under  given  operating  conditions,  these  polarization  losses  are  largely  dependent  on  cell  materials, 
electrode  microstructures,  and  cell  geometric  parameters.  Solid  oxide  fuel  cells  (SOFC)  with  yttria-stabilized  zirconia  (YSZ)  electrolyte, 
Ni-YSZ  anode  support,  Ni-YSZ  anode  interlayer,  strontium  doped  lanthanum  manganate  (LSM)-YSZ  cathode  interlayer,  and  LSM  current 
collector,  were  fabricated.  The  effect  of  various  parameters  on  cell  performance  was  evaluated.  The  parameters  investigated  were:  (1)  YSZ 
electrolyte  thickness,  (2)  cathode  interlayer  thickness,  (3)  anode  support  thickness,  and  (4)  anode  support  porosity.  Cells  were  tested  over  a 
range  of  temperatures  between  600  and  800  °C  with  hydrogen  as  fuel,  and  air  as  oxidant.  Ohmic  contribution  was  determined  using  the  current 
interruption  technique.  The  effect  of  these  cell  parameters  on  ohmic  polarization  and  on  cell  performance  was  experimentally  measured. 
Dependence  of  cell  performance  on  various  parameters  was  rationalized  on  the  basis  of  a  simple  analytical  model.  Based  on  the  results  of  the 
cell  parameter  study,  a  cell  with  optimized  parameters  was  fabricated  and  tested.  The  corresponding  maximum  power  density  at  800  °C  was 
~1,8W  cm~2. 
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1.  Introduction 

In  solid  oxide  fuel  cell  (SOFC)  research,  the  anode- 
supported  design  has  been  extensively  investigated  in  recent 
years.  This  is  because  anode-supported  SOFCs  can  be  fabri¬ 
cated  with  relative  ease,  are  mechanically  rugged,  and  exhibit 
some  of  the  highest  power  densities  at  temperatures  as  low 
as  750  °C.  Power  densities  in  excess  of  1  W  cm-2  at  the  op¬ 
erating  temperatures  have  been  reported  [1-3]. 

Although  significant  strides  have  been  made  in  recent 
years  in  improving  cell  performance,  there  is  still  consider¬ 
able  room  for  further  improvements,  even  with  the  standard 
materials  set  comprising  Y203-stabilized  zirconia  (YSZ)  as 
the  electrolyte  and  a  constituent  in  the  composite  anode  and 
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cathode,  Ni  for  the  anode,  and  Sr-doped  LaMnCE  (LSM)  for 
the  cathode.  Even  greater  improvements  are  deemed  possi¬ 
ble  with  other  materials  for  electrode  and  electrolyte  with 
intrinsically  superior  properties.  For  this  reason,  much  work 
is  currently  underway  in  this  area. 

There  are  five  loss  mechanisms,  which  result  in  the  loss  of 
useful  voltage  of  a  cell  at  a  given  operating  current  density. 
They  are:  (1)  ohmic  polarization,  which  is  the  voltage  loss 
due  to  the  ohmic  resistances  of  the  electrolyte,  electrodes  and 
interfaces  between  the  electrodes  and  the  electrolyte,  and  be¬ 
tween  the  electrodes  and  the  current  collectors  (contact  resis¬ 
tance).  (2)  Concentration  polarization  at  the  cathode,  which 
is  the  voltage  loss  associated  with  the  transport  of  gaseous 
oxidant  through  the  porous  cathode.  (3)  Activation  polariza¬ 
tion  at  the  cathode,  which  is  the  voltage  loss  associated  with 
the  oxygen  reduction  reaction.  (4)  Concentration  polariza¬ 
tion  at  the  anode,  which  is  the  voltage  loss  associated  with 
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the  transport  of  gaseous  fuel  through  the  anode.  (5)  Activation 
polarization  at  the  anode,  which  is  the  voltage  loss  associated 
with  the  hydrogen  oxidation  reaction.  The  total  polarization, 
jjtotah  can  be  expressed  as 

*7  total  =  ^electrode  "T  L'lnnic 

—  ^/a.act  +  ha.con  +  he, act  "T  he, con  T"  hohrnic  (1) 

where  heiectrode  is  the  sum  of  activation  and  concentration 
polarizations  at  the  two  electrodes,  ha.act  and  ha.con  are  re¬ 
spectively,  activation  and  concentration  polarizations  at  the 
anode,  hc.act  and  hc.con  are  respectively,  activation  and  con¬ 
centration  polarizations  at  the  cathode,  and  hohmic  is  the  total 
ohmic  polarization.  An  efficient  operation  of  SOFC  requires 
that  all  of  these  losses  be  as  small  as  possible.  Parameters  such 
as  the  ionic  conductivity  of  the  electrolyte  and  electronic  and 
ionic  conductivities  of  the  two  electrodes,  the  thicknesses  of 
the  electrolyte  and  electrodes,  and  possible  ohmic  resistances 
associated  with  interfaces,  determine  the  ohmic  loss.  It  is  of¬ 
ten  assumed  that  most  of  the  ohmic  loss  is  due  to  the  elec¬ 
trolyte.  Thus,  a  high  ionic  conductivity  and  a  small  electrolyte 
thickness  are  the  desired  characteristics  of  the  solid  elec¬ 
trolyte  to  minimize  the  ohmic  contribution.  Although  various 
solid  electrolytes  with  high  ionic  conductivities  at  moderate 
temperatures  have  been  explored,  yttria-stabilized  zirconia 
(YSZ)  is  by  far  the  most  widely  used  solid  electrolyte  due  to 
its  excellent  stability  in  both  reducing  and  oxidizing  environ¬ 
ments,  even  though  its  conductivity  is  lower  than  other  mate¬ 
rials  such  as  ceria  and  Sr-  and  Mg-doped  LaGaOs  (LSGM). 
Higher  conductivity  electrodes  and  lower  sheet  resistance  are 
also  important  in  lowering  the  ohmic  contribution.  Insofar 
as  electrodes  are  concerned,  the  ohmic  contribution  also  de¬ 
pends  upon  the  relative  amounts  of  the  two  phases  present,  the 
amount  of  porosity,  and  microstructure — in  particular  the  ge¬ 
ometry  of  the  particle  to  particle  contact.  Parameters  such  as 
the  electrode  porosity,  the  pore  size  and  the  pore  morphology 
(which  affects  the  totuosity  factor)  influence  the  transport  of 
gaseous  species  through  the  electrodes  and  thus  the  accompa¬ 
nying  concentration  polarization.  Finally,  the  electrode  inter¬ 
layer  morphology,  including  the  three  phase  boundary  (TPB) 
length,  determine  the  activation  polarization. 

In  addition  to  these  electrode  microstructure-related  pa¬ 
rameters,  the  electrode  thickness  is  also  expected  to  affect 
concentration  polarization;  the  thicker  is  the  support  elec¬ 
trode,  the  greater  is  the  concentration  polarization.  From  the 
standpoint  of  ruggedness  and  mechanical  integrity,  the  elec¬ 
trode  should  be  of  a  sufficient  thickness.  From  the  stand¬ 
point  of  concentration  polarization,  however,  the  thinner  is 
the  electrode,  the  lower  is  the  concentration  polarization.  It 
is  thus  necessary  to  determine  how  thick  the  electrode  sup¬ 
port  should  be  without  significantly  increasing  concentration 
polarization  losses. 

Within  a  given  set  of  materials,  such  as  YSZ  electrolyte, 
Sr-doped  LaMnOs  (LSM)  cathode,  and  Ni  +  YSZ  anode, 
large  variability  in  performance  characteristics  can  occur  de¬ 
pending  upon  the  geometric  design  of  the  fuel  cell  and  the 


microstructural  and  compositional  parameters  of  the  elec¬ 
trodes.  For  the  present  study  the  materials  selected  were 
LSM  of  composition  Lao.85Sro.i5MnO(3_s)  for  the  cath¬ 
ode  interlayer  and  cathode  current  collector,  Ni  for  the  an¬ 
ode  interlayer  and  anode  support,  and  YSZ  of  composition 
8  mol.%  Y2O3-92  mol.%  Zr02,  for  the  electrolyte,  for  the 
LSM  +  YSZ  cathode  interlayer,  for  the  Ni  +  YSZ  anode  in¬ 
terlayer,  and  for  the  Ni  +  YSZ  anode  support.  The  cathode 
interlayer,  where  most  of  the  electrocatalysis  related  to  oxy¬ 
gen  reduction  occurs  was  a  porous  composite  of  LSM  +  YSZ. 
The  anode  interlayer,  where  most  of  the  electrocatalysis  re¬ 
lated  to  hydrogen  oxidation  occurs  was  a  porous  composite 
of  Ni  +  YSZ.  The  relative  proportions  of  LSM  and  YSZ  in 
cathode  interlayer,  and  Ni  and  YSZ  in  anode  interlayer  and 
anode  support,  were  fixed  in  all  experiments.  The  principal 
objective  of  the  present  work  was  to  experimentally  study  the 
effect  of  electrolyte  thickness,  anode  support  thickness,  an¬ 
ode  support  porosity,  and  cathode  interlayer  thickness  on  cell 
performance.  Button  cells  were  fabricated  with  one  param¬ 
eter  varied  at  a  time,  keeping  other  parameters  fixed.  Cells 
were  subsequently  tested  with  humidified  hydrogen  as  fuel 
and  air  as  oxidant. 


2.  Experimental  procedure 

2.1.  Cell  fabrication 

A  typical  cell  used  in  the  present  work  consisted  of  the  fol¬ 
lowing  five  distinct  layers,  (a)  Porous  Ni  +  YSZ  anode  sup¬ 
port.  (b)  Porous  Ni  +  YSZ  anode  interlayer,  (c)  Dense  YSZ 
electrolyte,  (d)  Porous  LSM  +  YSZ  cathode  interlayer,  (e) 
Porous  LSM  cathode  current  collector.  Cell  fabrication  pro¬ 
cedure  is  briefly  described  in  what  follows.  NiO  and  YSZ 
powders  obtained  from  commercial  vendors  were  mixed  in 
the  desired  ratio,  70  wt.%  NiO+  30  wt.%  YSZ,  to  which  car¬ 
bon  powder  was  added.  The  amount  of  carbon  added  was 
based  on  the  desired  porosity.  Pellets  of  approximately  3  cm 
diameter  were  die-pressed.  A  thin  layer  of  NiO  +  Y SZ,  which 
is  the  anode  interlayer,  was  drop-coated  on  one  surface  of 
the  disc.  After  drying,  a  thin  layer  of  YSZ  electrolyte  was 
drop-coated  on  the  same  surface.  The  resulting  pellets  were 
sintered  in  air  at  1400  °C  for  2  h.  This  procedure  leads  to 
the  formation  of  a  fully  dense  YSZ  electrolyte  layer,  sup¬ 
ported  on  NiO  +  YSZ  anode  interlayer/NiO  +  YSZ  support. 
After  sintering  the  disc  is  about  2.6  cm  in  diameter.  A  thin 
layer  of  LSM  +  YSZ  containing  ~50  wt.%  LSM  +  ~50  wt.% 
YSZ  cathode  interlayer  was  painted  on  the  YSZ  layer.  The 
cell  was  then  fired  at  1200  °C  for  2  h.  After  firing  a  layer  of 
LSM  cathode  current  collector  was  applied.  The  area  of  the 
cathode  was  2  cm2,  which  was  used  as  the  basis  for  the  current 
density  calculation.  The  cell  was  then  heated  to  1 100  °C  to  en¬ 
sure  that  the  cathode  was  well  formed,  while  still  maintaining 
a  significant  level  of  porosity.  This  completes  cell  fabrication. 
One  parameter  was  varied  at  a  time  to  study  the  effect  of  a 
given  parameter.  The  parameters  varied  were:  (1)  electrolyte 
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thickness,  varied  between  ~4  and  ~20  |Jim.  (2)  Anode  sup¬ 
port  thickness,  varied  between  0.5  and  2.45  mm.  (3)  Anode 
support  porosity,  varied  between  32  and  76%  and  (4)  cathode 
interlayer  thickness,  varied  between  6  and  105  pm.  Param¬ 
eters  of  standard  cells,  against  which  performance  of  other 
cells  was  compared,  were  as  follows:  YSZ  electrolyte  thick¬ 
ness  pm,  anode  support  thickness  ~  1  mm,  anode  sup¬ 
port  porosity  ~48%,  cathode  interlayer  thickness  ~20  pm, 
anode  interlayer  thickness  ~20  pm,  and  cathode  current  col¬ 
lector  thickness  ~50  pm.  In  the  present  work,  anode  inter¬ 
layer  thickness  and  cathode  current  collector  thickness  were 
not  varied.  Also,  the  relative  proportions  of  constituents  in 
a  given  layer  were  not  varied,  and  neither  were  microstruc¬ 
tures  and  porosities  (with  the  exception  of  the  anode  support, 
whose  porosity  was  varied  over  a  range  as  stated  above). 

2.2.  Cell  testing 

For  testing,  a  typical  cell  was  mounted  in  a  test  fixture, 
a  schematic  of  which  is  shown  in  Fig.  1.  Silver  wires  were 
connected  to  silver  meshes  pressed  against  the  LSM  cathode 
current  collector  and  the  Ni  +  YSZ  anode  support,  respec¬ 
tively.  Two  sets  of  silver  leads,  one  on  each  side,  were  used 
for  current  collection  and  voltage  measurement.  Care  was 
taken  to  ensure  that  on  a  given  side,  the  two  silver  leads  were 
welded  on  the  mesh  at  the  same  point.  The  cell  was  heated  to 
800  °C  while  circulating  a  mixture  of  10%  FP  +90%  N2  on 
the  anode  side,  and  air  on  the  cathode  side.  NiO  from  the  an¬ 
ode  was  reduced  to  Ni  in  a  few  minutes  to  an  hour,  leading  to 
the  formation  of  additional  porosity.  Subsequently,  humid¬ 
ified  hydrogen  was  circulated  past  the  anode,  and  air  was 
circulated  past  the  cathode.  The  hydrogen  and  air  flow  rates 
were  300  and  550  ml  min- 1 ,  respectively.  Voltage  versus  cur¬ 
rent  density  polarization  curves  were  obtained  over  a  range 
of  temperatures  between  600  and  800  °C. 


Fig.  2.  An  oscilloscope  trace  illustrating  the  time  dependence  of  the  electric 
potential  before  and  after  current  interruption. 


2.3.  Measurement  of  ohmic  contribution 

Current  interruption  technique  was  used  to  obtain  the 
ohmic  contribution.  The  equipment  used  includes  a  Solartron 
SI  1287  Electrochemical  Interface  and  Agilent  54622 A  Dig¬ 
ital  Oscilloscope.  Fig.  2  shows  a  typical  voltage  versus  time 
trace  immediately  following  current  interruption.  The  sharp 
change  in  voltage  corresponds  to  the  ohmic  loss,  and  the  slow 
change  corresponds  to  non-ohmic  polarization  losses  (activa¬ 
tion  and  concentration).  Before  the  current  was  interrupted, 
care  was  taken  to  ensure  the  attainment  of  steady  state.  In  a 
few  cases,  current  interruption  tests  were  conducted  by  vary¬ 
ing  the  applied  current,  which  was  varied  between  0.5  and 
2  A.  Fig.  3  shows  the  measured  voltage  drop  as  a  function 
of  imposed  current  density.  A  linear  relation  between  volt¬ 
age  change  and  current  suggests  that  the  ohmic  portion  was 
accurately  captured  in  the  experiments.  Alternatively,  these 
experiments  suggest  that  the  time  constants  for  non-ohmic 
polarizations  were  considerably  longer  than  the  timescale  of 
the  current  interruption  experiment. 


Fig.  1 .  A  schematic  of  the  single  cell  testing  apparatus  used. 
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only  a  few  isolated  pores.  The  porous  regions  of  the  cathode 
interlayer,  the  anode  interlayer  and  the  anode  support  struc¬ 
ture  are  clearly  visible  in  the  figure.  Although  not  visible  in 
the  micrograph,  separate  measurements  have  shown  that  on 
the  anode  side,  Ni,  YSZ  and  porosity  form  contiguous,  three- 
dimensional  interpenetrating  networks,  and  there  is  a  consid¬ 
erable  amount  of  Ni-YSZ-pore  three-phase  boundary  (TPB) 
present.  The  existence  of  TPB  is  critical,  as  these  are  the  phys¬ 
ical  locations  where  the  electrochemical  reaction  of  hydrogen 
oxidation  occurs.  Similar  interpenetrating  networks  exist  be¬ 
tween  LSM-YSZ-pores  with  TPB  in  the  cathode  interlayer, 
where  the  electrochemical  reaction  of  oxygen  reduction  oc¬ 
curs.  Recent  work  has  shown  a  profound  effect  of  TPB  length 
in  composite  cathodes  on  cell  performance  [5]. 


3.2.  Single  cell  test  results 


Fig.  3.  Measured  ohmic  voltage  loss  as  a  function  of  current  density  on  a 
standard  cell  at  800  °C. 

2.4.  Microstructure  characterization 

After  testing,  some  of  the  cells  were  fractured  and  cross 
sections  were  examined  under  a  scanning  electron  micro¬ 
scope  (SEM).  The  fractured  pieces  of  some  of  the  samples 
were  vacuum-impregnated  with  an  epoxy.  After  hardening 
the  epoxy,  the  samples  were  polished  down  to  a  1  |Jim  finish. 
The  porosity  in  the  anode  interlayer,  anode  support,  cathode 
interlayer  and  cathode  current  collector  layer  was  determined 
by  quantitative  stereology  using  the  systematic  point  count 
procedure  [4], 

3.  Results 

3.1.  Microstructure  of  single  cells 

Fig.  4  shows  an  SEM  micrograph  of  the  fracture  surface 
of  a  typical  cell  after  testing.  Note  that  the  electrolyte  layer 
is  essentially  fully  dense,  as  evidenced  by  the  existence  of 


Fig.  4.  An  SEM  micrograph  of  a  typical  cell. 


Fig.  5  shows  the  performance  of  a  standard  cell  tested  at 
600,  700  and  800  °C.  The  maximum  power  density  (MPD) 
is  —0.2  W  cm-2  at  600  °C,  ~0.6  W  cm”2  at  700  °C  and 
—  1.2  W  cm-2  at  800  °C.  Fig.  6  shows  the  measured  volt¬ 
age  versus  current  density  traces  at  800  °C  for  cells  in  which 
the  only  parameter  varied  was  the  YSZ  electrolyte  thickness, 
which  was  varied  between  4  and  20  p,m,  all  other  parameters 
having  been  fixed  at  their  standard  values.  As  shown  in  Fig.  6, 
cells  with  8  and  15  |xm  electrolyte  thickness  exhibit  MPD  of 
about  1.2  W  cm-2,  while  that  with  20  |xm  exhibits  an  MPD 
of  1  W  cm-2.  The  cell  with  a  4  p,m  electrolyte  thickness  ex¬ 
hibits  MPD  of  —1.1  W  cm“2.  Fig.  7  shows  the  performance 
curves  for  the  same  cells  at  700  °C.  The  cell  with  an  8  |xm 
electrolyte  thickness  exhibits  MPD  of  —0.6  W  cm-2,  which 
is  much  lower  than  at  800  °C.  This  is  due  to  the  low  ionic 
conductivity  of  electrolyte  material  and  lower  rate  of  electro¬ 
chemical  reactions  at  lower  temperatures. 


Fig.  5.  Voltage  and  power  density  vs.  current  density  plots  for  a  standard 
cell  at  600,  700  and  800  °C. 
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Fig.  6.  Voltage  and  power  density  vs.  current  density  plots  at  800  °C  for 
cells  with  different  electrolyte  thicknesses. 


Fig.  8.  Voltage  and  power  density  vs.  current  density  plots  at  800  °C  for 
cells  with  different  anode  thicknesses. 


Fig.  8  shows  the  results  of  cell  tests  at  800  °C  wherein  the 
anode  support  thickness  was  varied  between  0.5  and  2.45  mm, 
all  other  parameters  fixed  at  the  standard  values.  Note  that  the 
MPD  varies  between  ~0.7  W  cm-2  for  anode  support  thick¬ 
ness  of  2.45  mm  to  ~1.35  W  cm~2  for  anode  support  thick¬ 
ness  of  0.5  mm  showing  a  profound  effect  of  anode  support 
thickness  on  performance.  Fig.  9  shows  the  results  of  cell  tests 
wherein  the  anode  support  porosity  was  varied  between  32 
and  76%,  all  other  parameters  maintained  at  the  standard  val¬ 
ues.  Note  that  the  MPD  for  a  cell  with  anode  support  porosity 
of  32%  is  only  ~0.72  W  cm-2,  while  that  for  a  cell  with  an¬ 
ode  support  porosity  of  57%  is  ~1.55  W  cm-2.  This  shows 
the  profound  effect  of  anode  support  porosity  on  cell  per¬ 
formance.  The  MPD  for  the  cell  with  anode  support  poros¬ 
ity  of  76%  was  about  1.5  W  cm-2.  The  lower  performance 
of  this  cell  despite  higher  porosity  of  the  anode  support  is 
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Fig.  7.  Voltage  and  power  density  vs.  current  density  plots  at  700  °C  for 
cells  with  different  electrolyte  thicknesses. 


due  to  the  fact  that  the  open  circuit  voltage  (OCV)  of  this 
cell  was  lower,  indicating  that  there  was  presumably  leakage 
through  the  electrolyte.  Fig.  10  shows  the  results  of  cell  tests 
with  all  parameters  maintained  at  the  standard  values,  with 
only  the  cathode  interlayer  thickness  varied  between  ~6  and 
~105  pirn.  Note  that  the  MPD  is  the  highest  for  cathode  inter¬ 
layer  thickness  of  20  pun,  and  is  lower  for  cells  with  cathode 
interlayer  thicknesses  either  smaller  or  larger  than  20  pim. 

Fig.  11  shows  the  results  of  cell  tests  with  optimized 
parameters:  anode  support  thickness  =  0.5  mm,  anode  in¬ 
terlayer  thickness  =  20  pim,  anode  support  porosity  =  57%, 
electrolyte  thickness  =  8  pirn,  cathode  interlayer  =  20  pim  and 
cathode  current  collector  =  50  pim.  In  the  present  work,  an¬ 
ode  interlayer  thickness  and  cathode  current  collector  thick¬ 
ness  were  not  varied.  The  MPD  is  ~1.8Wcm“2  at  800°C, 
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Fig.  10.  Voltage  and  power  density  vs.  current  density  for  cells  with  cathode 
interlayer  thickness  varied  between  ~6  and  ~  105  |xm. 


Fig.  12.  Measured  ohmic  area  specific  resistance  (ASR)  as  a  function  of 
YSZ  electrolyte  thickness  at  800  °C. 


~0.8  W  cm-2  at  700  °C,  and  -0.4  W  cm"2  at  600  °C.  The  re¬ 
sults  demonstrate  that  the  cell  with  the  optimized  parameters 
exhibits  the  highest  performance. 

3.3.  Area  specific  ohmic  resistance  (ASR) 

The  results  of  area  specific  ohmic  resistance  (ASR)  mea¬ 
sured  by  current  interruption  on  cells  as  a  function  of  elec¬ 
trolyte  thickness,  anode  support  thickness,  anode  support 
porosity  and  cathode  interlayer  thickness  are  given  in  Table  1 . 
As  seen  in  Table  1,  note  that  the  ASR  increases  with  increas¬ 
ing  electrolyte  thickness,  increasing  anode  support  thick¬ 
ness,  and  increasing  cathode  interlayer  thickness.  However, 
the  ASR  decreases  with  increasing  anode  support  porosity. 
Fig.  12  is  a  plot  of  ASR  at  800  °C  as  a  function  of  electrolyte 
thickness.  The  slope  of  the  plot  is  —24  12  cm  and  the  intercept 


Fig.  11.  Voltage  and  power  density  vs.  current  density  for  an  optimized  cell 
over  a  range  of  temperatures  between  600  and  800  °C. 


is  —0.087  £2  cm2.  Fig.  13  is  a  similar  plot  for  cell  tests  con¬ 
ducted  at  700  °C,  with  slope  equal  to  —63  12  cm  and  intercept 
—0. 138  ^2  cm2.  Fig.  14  is  a  plot  of  ASR  at  800  °C  as  a  function 
of  anode  support  thickness.  The  corresponding  slope  and  in¬ 
tercept  are  respectively,  0.24  12  cm  and  0.084  12  cm2.  Fig.  15 
is  a  plot  of  ASR  at  800  °C  as  a  function  of  cathode  interlayer 
thickness.  The  corresponding  slope  and  intercept  are  respec¬ 
tively,  3.922  12  cm  and  0.095  £2  cm2.  Finally,  Fig.  16  is  a  plot 
of  the  measured  ASR  at  800  °C  as  a  function  of  anode  support 
porosity. 

The  ASR  measured  by  current  interruption  on  standard 
cells  as  a  function  of  temperature  is  given  in  Table  2.  The 
ASR  increases  from  0.104  to  0.19  ( cm2  with  the  decrease  of 
temperature  from  800  to  700  °C.  The  ASR  for  the  cell  with 
optimized  parameters  is  also  given  in  Table  2.  Note  that  for 
the  optimized  cell,  the  ASR  increases  from  0.085  12  cm2  at 
800  °C  to  0. 15  12  cm2  at  700  °C. 
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Table  1 

Ohmic  resistance  of  cells  measured  by  current  interruption 

Anode  support  thickness  1  mm,  anode  interlayer  thickness  20  |xm,  electrolyte  thickness  8  |xm,  anode  support  porosity  48%,  cathode  interlayer  thickness 
20  |xm,  cathode  current  collector  thickness  50  fxm 

Electrolyte  thickness  (|xm)  4  8  15  20 

Ohmic  resistance  (ohm  cm2)  0.1  0.104  0.114  0.14 

Ohmic  resistance3  (ohm  cm2)  0.17  0.19  0.27  0.28 

Anode  support  porosity  48%,  anode  interlayer  thickness  20  [xm,  electrolyte  thickness  8  |xm,  cathode  interlayer  thickness  20  pm,  cathode  current  collector 
thickness  50  |xm 

Anode  support  thickness  (mm)  0.5  1.0  1.5  2.45 

Ohmic  resistance  (ohm  cm2)  0.095  0.104  0.13  0.14 

Anode  support  thickness  1  mm,  anode  interlayer  thickness  20  |xm,  electrolyte  thickness  8  |xm,  cathode  interlayer  thickness  20  |xm,  cathode  current  collector 
thickness  50  |xm 

Anode  support  porosity  (%)  32  48  57  76 

Ohmic  resistance  (ohm  cm2)  0.156  0.104  0.09  0.074 

Anode  support  thickness  1  mm,  anode  support  porosity  48%,  anode  interlayer  thickness  20  pm,  electrolyte  thickness  8  |xm,  cathode  current  collector 
thickness  50  |xm 

Cathode  interlayer  thickness  (|xm)  6  20  56  105 

Ohmic  resistance  (ohm  cm2)  0.095  0.104  0.12  0.135 

a  The  ohmic  resistance  measured  at  700  °C.  All  other  measurements  at  800  °C. 
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Fig.  14.  Measured  ohmic  area  specific  resistance  (ASR)  as  a  function  of 
anode  support  thickness  at  800  °C. 
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Fig.  15.  Measured  ohmic  area  specific  resistance  (ASR)  as  a  function  of 
cathode  interlayer  thickness  at  800  °C. 


Fig.  16.  Measured  ohmic  area  specific  resistance  (ASR)  as  a  function  of 
anode  support  porosity  at  800  °C. 


Table  2 


Ohmic  resistance  (by  current  interruption) 
standard  and  optimized  cells 

as  a  function  of  temperature  for 

Temperature  (°C) 

700 

800 

Ohmic  ASR  of  standard  cell  (ohm  cm2) 

0.19 

0.104 

Ohmic  ASR  of  optimized  cell  (ohm  cm2) 

0.15 

0.085 

3.4.  Measurement  of  porosity 

Measured  porosity  levels  using  quantitative  stereology 
were  ~26%  for  the  cathode  interlayer,  ~45%  for  the  cathode 
current  collector,  and  ~23%  for  the  anode  interlayer. 
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4.  Discussion 


4.1.  Effect  of  various  parameters  on  the  ohmic  area 
specific  resistance  (ASR) 


For  a  given  set  of  materials,  compositions  and  microstruc¬ 
tures,  the  ASR  is  a  function  of  the  thicknesses  of  the  elec¬ 
trolyte,  the  two  electrodes,  and  the  possible  interfacial  or  con¬ 
tact  resistances.  In  the  present  work,  each  cell  comprised  of 
five  distinct  layers;  porous  Ni  +  YSZ  anode  support,  porous 
Ni  +  YSZ  anode  interlayer,  dense  YSZ  electrolyte,  porous 
LSM  +  YSZ  cathode  interlayer,  and  porous  LSM  current  col¬ 
lector.  Thus,  the  ASR,  //,.  may  be  given  as 


n  ionic  i  ,  elect  /  ,  elect  j  ,  elect; 

Ri  —  Pe  ie  +  Pc(1)  Zc(l)  +  Pc(2)  4(2 )  +  Pa( p  4(1) 

T  f>,,n  i  l at 2)  "b  ^contact 


where  p",nlc  and  Ze  are  respectively,  the  electrolyte  ionic  re¬ 
sistivity  and  thickness,  p®^1  and  ZC( i)  are  respectively,  the 
cathode  current  collector  electronic  resistivity  and  thickness, 
b’c'zf  an<^  4(2)  are  respectively,  cathode  interlayer  electronic 
resistivity  and  thickness,  paJ“l  and  Za(i)  are  respectively,  an¬ 
ode  support  electronic  resistivity  and  thickness,  pf^f  and  4(2) 
are  respectively,  anode  interlayer  electronic  resistivity  and 
thickness,  and  RContact  is  some  contact  resistance,  representa¬ 
tive  of  resistance  associated  with  interfaces.  A  schematic  of  a 
typical  cell  with  various  regions  labeled  is  shown  in  Fig.  17. 

In  cases  wherein  the  only  parameter  varied  is  the  thickness 
of  one  of  the  layers,  the  thickness  dependence  of  resistance 
can  be  unequivocally  attributed  to  the  resistivity  of  that  layer. 
For  example,  in  cells  in  which  the  only  parameter  varied  was 
the  YSZ  electrolyte  thickness,  the  dependence  of  ASR  on 
thickness  can  be  uniquely  related  to  the  resistivity  of  dense 
YSZ.  It  is  expected  that  cells  with  the  thinnest  electrolyte 
(~4  p,m)  should  exhibit  the  highest  performance.  The  obser¬ 
vation  that  the  cell  with  the  ~4  p,m  electrolyte  does  not  ex¬ 
hibit  the  highest  MPD  as  shown  in  Fig.  6  simply  implies  that: 


(1)  the  cells  were  not  identical  in  all  other  respects,  despite 
the  intent.  For  example,  there  might  have  been  small  (unin¬ 
tended)  differences  in  the  thicknesses  and/or  microstructures 
of  the  other  layers.  (2)  The  electrolyte  ohmic  contribution  is 
not  the  overriding  factor  at  these  low  electrolyte  thicknesses. 

Since  the  only  parameter  varied  is  the  electrolyte  thick¬ 
ness,  the  Eq.  (2)  can  be  written  as 

R  i  =  Pr%  +  R  const  (3) 

where  //const  is  the  contribution  to  ASR  from  sources  other 
than  the  electrolyte.  Thus,  a  plot  of  the  measured  ASR  versus 
the  electrolyte  thickness,  Ze,  should  be  linear  with  the  slope 
equal  to  the  electrolyte  ionic  resistivity,  p"’"lc  and  intercept 
Rconst.  A  plot  of  the  ohmic  area  specific  resistance  (ASR)  at 
800  °C  measured  by  current  interruption  versus  Ze  is  shown 
in  Fig.  12  with  the  slope  ~24  Q  cm,  which  is  in  good  agree¬ 
ment  with  the  reported  value  of  the  resistivity  of  YSZ  elec¬ 
trolyte  (of  typical  grain  size  on  the  order  of  a  few  microns)  of 
22.3  £2  cm  at  800  °C  [6].  Note,  however,  that  the  intercept  is 
nonzero  (0.087  cm2),  indicating  that  there  are  other  sources 

of  substantial  ohmic  contribution.  Fig.  13  shows  a  similar 
plot  of  the  measured  ASR  by  current  interruption  versus  Ze  at 
700  °C.  The  measured  electrolyte  resistivity  from  the  plot  is 
~63  Q  cm  and  the  corresponding  intercept  is  ~0. 1 38  Q  cm2, 
showing  that  a  substantial  contribution  to  the  ohmic  contri¬ 
bution  from  sources  other  than  the  electrolyte  is  also  present 
at  700  °C. 

Fig.  14  is  a  plot  of  the  ASR  as  a  function  of  anode  support 
thickness  at  800  °C.  Since  the  only  parameter  varied  was  the 
anode  support  thickness,  the  ASR  as  a  function  of  anode 
support  thickness  is  given  by 

—  Pa(l)l4(l)  +  ^const  (4) 

where  the  slope  is  the  electronic  resistivity  of  the  anode  sup¬ 
port,  and  the  intercept  is  Rc0nst,  which  is  the  ohmic  part 
due  to  sources  other  than  the  anode  support.  From  Fig.  14, 
the  estimated  value  of  paj“l  is  ~0.24  cm  and  the  intercept, 


Fig.  17.  A  schematic  diagram  of  an  anode-supported  cell  comprising  five  distinct  layers,  and  variation  of  partial  pressures  of  the  various  gaseous  species. 
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Table  3 

Calculated  electrolyte  ohmic  contribution  and  measured  total  ohmic  contribution  (by  current  interruption)  at  800  °C 


YSZ  thickness  (|xm) 

^electrolyte  (calculated)2  (£2  cm2) 

^ohmic  (measured)  (£2  cm2) 

A^ohmic  (£2  Cm2) 

4 

0.0096 

0.10 

0.0904 

8 

0.0192 

0.104 

0.0848 

15 

0.036 

0.114 

0.078 

20 

0.048 

0.140 

0.092 

a  Electrolyte  resistivity  =  24  Q  cm. 


^const  's  ~0.084£2cm2.  In  the  anode  support,  the  volume 
fraction  of  Ni  is  ~35%,  the  remainder  being  porosity  (essen¬ 
tially  infinite  resistivity)  and  YSZ  (also  of  very  high  resistiv¬ 
ity,  in  comparison  to  Ni).  That  is,  for  all  practical  purposes, 
the  volume  fraction  of  conductive  phase  (Ni)  is  ~35%.  It  is 
known  that  at  relatively  high  volume  fractions  of  insulating 
phases,  the  net  conductivity  can  be  far  below  that  of  the  con¬ 
ductive  phase  corrected  for  porosity  through  a  simple  linear 
relationship,  provided  the  inter  particle  contact  is  poor.  For 
example,  it  has  been  found  that  a  plot  of  relative  conductiv¬ 
ity,  cr( V„)/er0,  versus  relative  density,  £(Vv)/£, 0  where  er0  is  the 
conductivity  of  a  fully  dense  material,  er(V„)  is  the  conduc¬ 
tivity  of  a  material  with  Vv  volume  fraction  porosity,  £0  is 
the  density  of  fully  dense  material,  and  t;(Vv )  is  the  density 
of  a  material  with  porosity  Vv  is  often  approximately  linear 
and  extrapolation  intersects  the  axis  at  ~0.45-~0.6. 

That  is,  the  er(Vv)/er0  almost  approaches  zero  for  Vv  >  0.45 
[7-9].  That  is,  often  the  contiguity  between  conductive  par¬ 
ticles  is  very  poor.  In  light  of  this,  even  though  er0  for  Ni  at 
800°C  is  very  high  (~7.94  x  104  S  cm-1),  it  is  not  surpris¬ 
ing  that  the  estimated  er( Vv)  (for  Vv  ~  0.65)  is  only  1/0.24  or 
~4.4S  cm-1. 

Fig.  15  is  a  plot  of  ASR  as  a  function  of  cathode  interlayer 
thickness  at  800  °C.  Since  the  only  parameter  varied  was  the 
cathode  interlayer  thickness,  the  ASR  as  a  function  of  cathode 
interlayer  thickness  is  given  by 

=  Pc(2)  ^c(2)  +  ^const  (5) 

where  the  slope  is  essentially  the  electronic  resistivity  of  the 
cathode  interlayer,  pJ:J“l,  and  the  intercept  is  R"onst,  which  is 
the  ohmic  contribution  from  sources  other  than  the  cathode 
interlayer.  As  seen  in  Fig.  15,  the  slope,  p^f  is  estimated 
to  be  3.922  Q  cm  and  the  intercept,  R"onst  is  ~0.095  £2  cm2. 
The  electronic  resistivity  of  fully  dense  LSM  at  800  °C  is 
~2  x  10-2  £2  cm.  The  fact  that  the  estimated  p®^)1  from  the 
slope  is  considerably  higher  (lower  conductivity)  is  consistent 
with  the  rather  high  volume  fraction  of  the  insulating  phase. 


namely  porosity  +  YSZ  (Vv~0.65),  such  that  connectivity 
(particle  to  particle  contact)  is  poor. 

Independent  measurement  of  cell  ASR  as  a  function  of 
relative  thicknesses  thus  has  facilitated  the  estimation  of  the 
various  resistivities,  namely,  Pg0nic,  p^1,  and  p^1.  In  the 
standard  cells,  the  le,  /a(  | ,  and  lcl2)  were  respectively,  8  |xm, 
1  mm,  and  20  |xm.  Thus,  for  the  standard  cells,  the  estimated 
ohmic  ASR  at  800  °C  is 

Ri  =  24  x  8  x  10-4  +  0.24  x  10“ 1  +  3.922  x  20  x  10-4 

+  Cnst  =  0-051 +  /CnstfW 

wherein 

-^const  =  Pa(2)^a(2)  +  Pc(l)^c(l)  +  ^contact 

The  experimentally  measured  ASR  by  current  interrup¬ 
tion  is  ~0. 1 04  T2  cm2.  Thus,  the  estimated  value  of  R"'nst 
is  0.104  —  0.051  =0.053  £2  cm2.  That  is,  approximately  half 
of  the  ohmic  ASR  at  800  °C  of  the  standard  cells  used  in 
this  work  is  attributed  to  resistances  of  the  Ni  +  YSZ  an¬ 
ode  interlayer,  LSM  cathode  current  collector,  and  any  pos¬ 
sible  contact  resistances.  Tables  3  and  4,  respectively,  list 
the  measured  cell  ohmic  ASR  at  800  and  700  °C,  along  with 
calculated  electrolyte  contribution  for  electrolyte  thickness 
ranging  between  4  and  20  p,m.  Note  that  even  for  an  elec¬ 
trolyte  thickness  of  20  |xm  and  at  700  °C,  its  ohmic  contri¬ 
bution  to  the  ASR  is  lower  than  that  due  to  other  sources. 
For  a  8  p,m  electrolyte  thickness,  the  electrolyte  contribution 
to  the  ohmic  ASR  at  800  °C  is  ~0.0192  £2  cm2,  which  cor¬ 
responds  to  ~19%  of  the  total  ohmic  ASR.  At  700  °C,  the 
electrolyte  contribution  is  ~0.0504  £2  cm2,  which  is  ~27% 
of  the  total  ohmic  ASR.  This  shows  that  even  at  temper¬ 
atures  as  low  as  700  °C,  the  ohmic  contribution  of  a  thin 
(~8  |xm)  YSZ  electrolyte  is  rather  small,  and  that  YSZ  is  a 
satisfactory  electrolyte.  Also,  it  is  seen  that  there  is  no  need 
to  lower  the  electrolyte  thickness  much  below  about  8  |xm, 
since  the  electrolyte  contribution  to  the  overall  ASR  is  rather 
modest.  That  is,  although  an  electrolyte  material  exhibiting 


Table  4 

Calculated  electrolyte  contribution  and  measured  total  ohmic  contribution  (by  current  interruption)  at  700  °C 


YSZ  thickness  (|xm) 

^electrolyte  (calculated)2  (£2  cm2) 

^ohmic  (measured)  (£2  cm2) 

A^ohmic  (£2  Cm2) 

4 

0.0252 

0.17 

0.145 

8 

0.0504 

0.19 

0.14 

15 

0.0945 

0.21 

0.12 

20 

0.126 

0.28 

0.154 

Electrolyte  resistivity  =  63  £2  cm. 
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higher  ionic  conductivity  than  YSZ  is  preferred,  it  is  not  an 
absolute  necessity  for  intermediate  temperature  SOFC  oper¬ 
ating  at  temperatures  >700  °C.  This  also  means  that  conven¬ 
tional,  thick  film  processing  methods  capable  of  producing 
supported  YSZ  membranes  of  a  thickness  on  the  order  of  a 
few  microns,  are  satisfactory  for  the  fabrication  of  electrolyte 
films  of  SOFC.  Further  work  will  be  required  to  determine 
^ajzf  anc*  Pc  i?  by  varying  respectively,  the  anode  interlayer 
thickness  and  the  cathode  current  collector  thickness.  If  this 
is  done,  it  will  allow  the  estimation  of  the  contact  resistance, 

^contact  ■ 

Fig.  16  shows  the  dependence  of  ohmic  ASR  as  a  function 
of  anode  support  porosity.  Over  the  range  of  porosities  inves¬ 
tigated,  the  ASR  decreases  with  increasing  porosity — from 
0.156  £2  cm2  at  32%  porosity  to  0.074  £2  cm2  at  76%  poros¬ 
ity.  This  behavior  appears  unreasonable,  since  it  would  be  ex¬ 
pected  with  increase  in  porosity  (decrease  in  volume  fraction 
of  Ni),  the  net  contribution  to  the  ohmic  ASR  should  increase. 
The  volume  fraction  of  Ni  in  the  sample  containing  32% 
porosity  was  ~0.47,  and  that  containing  76%  porosity  was 
~0.16.  The  resistivity  of  Ni  at  800  °C  is  1.26  x  10-5  £2  cm 
or  conductivity  of  7.94  x  104  S  cm-1.  If  a  simple  linear  re¬ 
lationship  can  be  assumed  between  conductivity  and  vol¬ 
ume  fraction  of  Ni,  namely  er(UNj)  ^  rr0  UN[,  the  expected 
contributions  to  ASR  for  a  1  mm  thick  anode  support  are 
2.68  x  10-6  £2  cm2  for  32%  porosity  and  7.87  x  10-6  £2  cm2 
for  76%  porosity.  Clearly,  these  are  rather  small  values 
compared  to  the  experimental  measurements,  which  are  re¬ 
spectively,  0.156  and  0.074  £2  cm2,  or  the  measured  ASRs 
are  greater  by  factors  of  0.156/(2.68  x  10_6)~  58,210  and 
0.074/(7.87  x  10_6)~9400),  respectively.  This  means  that 
there  must  be  other  factors  indirectly  dependent  on  poros¬ 
ity,  which  exhibit  greater  effect  on  the  overall  conductiv¬ 
ity  of  anode  support.  The  most  likely  one  is  the  nature  of 
contiguity — namely  the  inter  particle  contact  geometry.  An 
approximate  order  of  magnitude  effect  of  the  geometry  of 
contact  can  be  determined  by  making  simple  calculations 
based  on  an  assumed  geometry.  If  nickel  particles  are  as¬ 
sumed  to  be  spherical  of  radius  r  and  if  the  inter  particle  neck 
radius  is  r(l,  it  can  be  shown  that  the  effective  resistivity  of  a 
body  comprising  spheres  is  given  by 


Peff 


”  lnj1  +  ^zf:| 

2n/1  -  X1  l  1  -  Vl  -k2  j 


(6) 


where  X  =  r0lr ,  is  a  measure  of  the  relative  neck  size — or  al¬ 
ternatively,  how  good  the  contact  is  between  particles.  Note 
that  as  X  ->  0,  the  effective  resistivity  of  the  body  approaches 
infinity,  showing  the  profound  effect  of  contact  geometry. 
Presumably,  the  microstructure  in  the  porous  electrodes  must 
be  such  that  in  the  sample  with  76%  porosity,  even  though  the 
volume  fraction  of  Ni  is  low  (~0. 16),  the  contiguity  must  be 
greater  than  the  sample  with  larger  Ni  content.  Unfortunately, 
in  the  absence  of  detailed  quantitative  microstructural  char¬ 
acterization  of  the  anode,  this  question  remains  unanswered 
at  the  present  time. 


4.2.  Concentration  and  activation  polarizations: 
comparison  of  experimental  results  on  cell  tests  with  an 
analytical  model 


Fig.  17  shows  a  schematic  of  the  cross  section  of  a  cell, 
along  with  schematic  variations  in  the  partial  pressures  of 
various  gaseous  species.  In  terms  of  the  various  parameters 
related  to  transport  and  electrode  reactions,  it  is  possible  to 
describe  the  voltage  versus  current  density  relationship  math¬ 
ematically.  Several  models  have  been  proposed,  with  vary¬ 
ing  degrees  of  complexity  and  requiring  a  number  of  differ¬ 
ent  parameters  [1,10,11].  In  this  work,  one  of  the  simplest 
models  requiring  a  relatively  small  number  of  parameters  is 
used.  The  partial  pressures  of  hydrogen  in  the  fuel  (outside 
of  the  anode)  and  of  oxygen  in  the  oxidant  (outside  of  the 
cathode)  are  given  respectively,  by  p^  and  p\'h .  The  fuel 
is  usually  humidified.  The  partial  pressure  of  H2O  in  the 
fuel,  just  outside  the  anode,  is  the  Present  work, 

the  values  of  hydrogen  and  FLO  partial  pressures  were  re¬ 
spectively,  p'f  0.97  atm  and  /?jji0  ~  0.03  atm.  The  partial 
pressures  of  hydrogen  and  FLO  at  the  anode  support/anode 
interlayer  interface  corresponding  to  a  current  density  of  i 
are  given  respectively,  by  pn^(i)O)  and  PH->O(i)0),  and  that 
of  oxygen  at  the  cathode  current  collector/cathode  interlayer 
interface  is  given  by  po2(i)(i)-  Finally,  partial  pressures  of  hy¬ 
drogen  and  H2O  at  the  anode  interlayer/electrolyte  interface 
are  given  respectively,  by  p'^ppO  and  PhiO(<)(*)’  anc*  that  °f 
oxygen  at  the  cathode  interlayer/electrolyte  interface  is  given 
by  ,.)(/).  In  what  follows,  it  is  assumed  that  the  flow  rates 
of  fuel  and  oxidant  are  sufficiently  fast  such  that  /?]  j  and  p'ft 
are  constant,  independent  of  current  density  (that  is,  both  fuel 
and  oxidant  utilization  are  negligible)  [12].  The  partial  pres¬ 
sures  of  hydrogen  and  oxygen  in  the  electrodes,  however,  are 
functions  of  current  density.  It  can  be  shown  that  that  anodic 
concentration  polarization  with  an  FL-FLO  gas  mixture  as 
fuel  is  given  by  [1] 


77  a,  con  — 


RT 

- In 

2  F 


/  /7h2(0^^',H2o\ 
\^H2^H20(i)(lV 


(7) 


where  R  is  the  ideal  gas  constant  and  F is  the  Faraday  constant. 
Local  equilibrium  is  assumed,  which  for  anode  means  that 
locally,  ph9,  P\\-,o,  and  /?q,  are  related  to  each  other  via  the 
equilibrium 

H2  +  \02  h2o 

In  what  follows,  it  will  be  assumed  that  gaseous  transport 
through  the  porous  electrodes  can  be  adequately  described 
by  binary  diffusion  and  that  there  is  negligible  contribution 
of  viscous  flow.  The  microstructure  within  a  given  region, 
e.g.  the  anode  support  or  the  anode  interlayer,  is  expected 
to  be  uniform.  Under  these  conditions,  the  partial  pressures 
of  FL  and  FLO  within  a  given  region  are  expected  to  vary 
linearly  with  position,  as  shown  in  Fig.  17.  It  is  to  be  noted, 
however,  that  the  partial  pressure  of  oxygen  at  the  anode, 
which  is  much  lower  than  the  partial  pressures  of  hydrogen 


F.  Zhao,  A.V.  Virkar  /  Journal  of  Power  Sources  141  (2005)  79-95 


89 


and  HtO,  i.e.  po2(a)  P\h  •  P\  h(h  does  not  vary  linearly,  but 
is  governed  by  the  above  Ho/bbO  equilibrium.  It  is  also  to  be 
mentioned  that  the  variation  of  p\\n  and  p\\-,o  in  the  anode 
interlayer,  especially  close  to  the  electrolyte,  is  not  expected 
to  be  linear.  This  is  because  within  this  layer,  individual  fluxes 
in  the  gas  phase  are  not  conserved,  even  though  p\\2  +  pu-,o 
continues  to  remain  fixed.  This  aspect  is  ignored  here. 

The  effective  anode  support  and  anode  interlayer  diffu- 
sivities  are  respectively,  and  wherein  it 

is  assumed  that  the  effective  diffusivities  are  proportional  to 
the  H2-H2O  binary  diffusivity  ( Dh2-h2o)>  volume  fraction 
porosity,  and  inversely  proportional  to  tortuosity  factor.  Since 
other  effects,  such  as  Knudsen  diffusion  are  likely  present, 
and  the  exact  nature  of  tortuosity  is  unclear,  in  the  above 
description  the  tortuosity  factor  is  a  merely  phenomenolog¬ 
ical  fitting  parameter.  In  steady  state,  under  the  assumptions 
made,  the  partial  pressures  of  hydrogen  and  H2O  at  the  an¬ 
ode  interlayer/electrolyte  interface  as  a  function  of  current 
density  are  given  respectively,  by 


Pn2(i)^ 
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H2-H20  L 
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and 
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The  maximum  possible  current  density  corresponds  to  the 
lowest  possible  partial  pressure  of  hydrogen  at  the  anode  in¬ 
terlayer/electrolyte  interface,  which  is  zero  (although  it  can 
never  be  exactly  zero).  Thus,  setting  Eq.  (8)  equal  to  zero 
gives  the  anode  limiting  current  density,  given  by 


eff(l) 


2  P[  h  \  It 
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,eff(2) 


H2  H2-H20^H2-H20 


7?  ?4(1)  4(2)  (  ^jj2 — h2  O  /  4(1 )  ' 
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If  all  of  the  requisite  parameters  are  known,  for  a  given 
current  density,  i,  the  and  Ph2o((')0)  can  he  estimated 

from  Eqs.  (8)  and  (9),  respectively,  and  the  anodic  concen¬ 
tration  polarization  can  be  calculated  using  Eq.  (7). 

Concentration  polarization  at  the  cathode  can  be  similarly 
examined.  The  oxidant  is  typically  air,  a  mixture  of  O2  and 
N2. 

A  schematic  of  a  porous  cathode,  along  with  the  spa¬ 
tial  variation  of  partial  pressure  of  oxygen  is  also  shown  in 
Fig.  17.  For  an  anode-supported  cell  with  a  small  cathode 
thickness,  the  oxygen  partial  pressure  at  the  interface  of  cur¬ 


rent  collector  and  cathode  interlayer  can  be  given  by  [13] 


Po2(i)(0 


P  02 


(p-P°o2\ 
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where  p°j  is  the  oxygen  partial  pressure  in  the  oxidant  (air) 
outside  of  the  cathode  (Pa  or  atm),  po2(i)(j)  is  the  oxygen 
partial  pressure  at  the  interface  between  the  cathode  current 
collector  and  the  cathode  interlayer  corresponding  to  cur¬ 
rent  density  i,  p  the  total  pressure,  4(1)  the  cathode  current 
collector  thickness  and  is  the  effective  O2-N2  diffu¬ 

sivity  through  the  cathode  current  collector.  The  oxygen  flux 
through  the  current  collector  and  cathode  interlayer  should 
be  the  same,  i.e. 


(Po2  -  Po2(i)(0)  I  ~~P~ 

\ 1  po2 

=  (po2(i)(i)  -  Po2(i)(0) 


neff(l) 
i)2  — No 

4(1) 

P 


P-PO2(i)(0j  4(2) 


D: 


,eff(2) 
02—  N2 


(12) 


where  Pq^Si)  is  the  oxygen  partial  pressure  at  the  interface 

between  the  cathode  interlayer  and  electrolyte,  Og^-No  the 
effective  O2-N2  diffusivity  through  the  cathode  interlayer, 
and  4(2)  is  the  cathode  interlayer  thickness.  Here  also,  the 
variation  of  po2  in  the  cathode  interlayer  will  be  nonlinear, 
especially  close  to  the  electrolyte/cathode  interlayer  inter¬ 
face,  as  the  flux  of  O2  is  not  conserved.  This  aspect  is  ignored 
here. 

Inserting  Eq.  (11)  into  Eq.  (12)  and  rearranging 
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The  maximum  possible  current  density  is  that  correspond¬ 
ing  to  the  lowest  possible  oxygen  partial  pressure  at  the  in¬ 
terface  between  the  cathode  interlayer  and  the  electrolyte, 
which  is  zero  (although  Pq^Ai)  of  course  can  never  be  ex¬ 
actly  zero).  The  corresponding  current  density,  which  is  the 
cathode  limiting  current  density,  ;cs,  is  given  by  setting  Eq. 
(13)  to  zero,  and  is  given  by 

-B+  y/B2  +4AC 


where  the  constants  A-C  are  given  by 
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and 
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The  concentration  polarization  at  the  cathode  is  given  by 
RT 


»7c,< 
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Finally,  the  voltage  versus  current  density,  namely,  V(i) 
versus  i  is  given  by  [1] 


RT 
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where  it  is  assumed  that  the  net  activation  polarization  at  both 
the  cathode  and  anode  can  be  described  by  the  Tafel  equation 
and  is  given  by 


Fig.  18.  Experimental  voltage  vs.  current  density  plot  at  800  °C  for  a  cell 
with  anode  support  porosity  of  32%  and  the  corresponding  best  fit  to  Eq. 
(19).  The  fit  is  shown  only  for  current  density  above  about  0.2  A  cm-2  since 
the  Tafel  equation  is  not  satisfactory  at  too  low  a  current  density. 


riact  =  a  +  b\ni  (20) 

Note  that  and  p'H20(i)(i)  are  functions  of 

and  and  p'0l(i)(i)  is  a  function  of  £>02-n2  and 

which  are  four  of  the  unknown  parameters.  The 
other  unknown  parameters  are  a  and  b,  which  describe  the 
overall  activation  polarization.  Strictly,  a  and  b  are  functions 
of  Ph,(;)(0  and  and  thus  themselves  are  functions 

of  the  current  density  i.  This  aspect  is  ignored  here  insofar  as 
curve-fitting  is  concerned.  However,  implications  of  their  de¬ 
pendence  on  partial  pressures  are  discussed  later.  The  ohmic 
ASR,  Rx,  is  known  from  the  current  interruption  tests.  The 
experimental  V(i)  versus  i  trace  for  each  case  was  fitted  us¬ 
ing  seven  adjustable  parameters:  Aj,  a ,  b, 

anc*  These  parameters  were  adjusted  to 

obtain  the  best  possible  fit  to  Eq.  (19).  In  the  curve-fitting 
procedure,  the  Aj  was  also  treated  as  an  adjustable  parameter. 
Fig.  18  shows  a  comparison  with  the  experimental  data  on 
a  cell  with  anode  support  porosity  of  32%  tested  at  800  °C 
and  the  best  fit  to  the  data.  Note  that  the  fit  is  very  good, 
except  at  very  low  current  densities,  where  the  fit  diverges. 
This  is  simply  the  result  of  assuming  the  applicability  of  the 
Tafel  equation  over  the  entire  range.  However,  the  Tafel  equa¬ 
tion,  which  has  a  term  with  In  i,  becomes  divergent  as  i  ->  0. 
Thus,  the  Tafel  equation  is  not  applicable  at  too  low  a  cur¬ 
rent  density  (lower  than  the  exchange  current  density,  i0),  and 
one  must  use  an  effective  charge  transfer  resistance  given  by 
Rct (eff)  ~  RT/zFia  with  the  activation  polarization  given  by 
?7act  ^  t  Rct(eff)  •  This  aspect  was  ignored  in  fitting  the  data; 
and  hence  the  divergence  of  the  fit  at  low  current  densities. 
That  part  of  the  fit,  therefore,  has  been  shown  by  a  dotted 
line.  Table  5  gives  the  parameters  A1,,  a,  b,  /}q!-n2’ 

anc*  corresponding  to  the  best  fits  to  all 

of  the  cell  tests  conducted. 


The  first  part  of  Table  5  gives  the  fitting  parameters  for 
cells  in  which  the  only  parameter  varied  was  the  anode  sup¬ 
port  thickness,  which  was  varied  between  0.5  and  2.45  mm, 
all  other  parameters  having  been  fixed  at  the  standard  val¬ 
ues.  The  second  column  gives  the  Rj,  which  varies  between 
0.101  and  0.148  £2  cm2.  Comparison  with  Table  1  shows  that 
the  ASR  determined  by  current  interruption  varies  between 
0.095  and  0. 14  £2  cm2,  in  good  agreement  with  the  Aj  ob¬ 
tained  through  fitting  to  the  model  (Eq.  (19)).  All  V(i)  versus 
i  traces  from  this  set  could  be  adequately  fitted  using  one 

.  f  i  f  r.eff(U  r>eff(2)  r>eff(0  A 
set  of  values  for  /\)2_NV  Oo2-n2’  «„2_h,0  and  A)n,_„20, 

which  were  respectively,  0. 14,  0.04,  0.68,  and  0.08  cm2  s“ 1 . 
Since  the  only  parameter  varied  was  the  anode  support 
thickness,  it  is  expected  that  the  fitted  transport  parame¬ 
ters  should  be  the  same,  consistent  with  the  fitting  results. 
Anode  support  and  cathode  current  collector  have  greater 
porosities  than  respectively,  anode  interlayer  and  cathode  in¬ 
terlayer.  This  is  consistent  with  >  ^nl-ibO  ar,d 

^Cb-No  >  an  earlier  study,  from  the  measure¬ 

ment  of  /as,  the  T^-h-iO  was  estimated  to  be  ~0.5  cm2/s 
[12].  In  that  study,  a  single  T>h!  [[,0  was  assumed  for  the 
entire  anode,  despite  the  presence  of  two  distinct  layers  (an¬ 
ode  interlayer  and  anode  support).  The  present  curve-fitting 
shows  that  the  effective  diffusivity  through  the  anode  inter¬ 
layer  (~0.08  cm2  s_1)  is  much  lower  than  through  the  anode 
support  (~0.68  cm2  s_  1 ),  mainly  due  to  its  much  lower  poros¬ 
ity  (~23%)  compared  to  that  in  the  anode  support  (~48%). 
The  present  estimates  of  ~0.68  cm2  s~ 1  for  and 

~0.08cm2s_1  for  are  dlus  reasonable.  Indepen¬ 

dent  measurements  using  an  electrochemical  technique  for 
LSM  porous  bodies  with  porosity  ranging  between  ~15  and 
~45%  showed  that  the  corresponding  ranged  be¬ 

tween  ~0.03  and  ~0.12cm2  s-1  [13].  The  estimated  values 
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Table  5 

Fitted  parameters  for  the  anode-supported  cells 


R[  (£2  cm2)  a  b 


io  (mA  cm  2)  0*%,  (cm2s  ') 


D 


etf(2) 

02-N2 


(cm2  s'* ) 


D 


etf(l) 

H2-H20 


(cm2  s'1) 


D 


etf(2) 

h2-h2o 


(cm2  s  1 ) 


Anode  support  thickness  varied  (electrolyte  thickness  =  8  |xm,  anode  support  porosity  =  48%,  cathode  interlayer  thickness  =  20  p,m) 
Anode  support  thickness  (mm) 


0.5 

0.101 

0.098 

0.11 

410 

0.14 

0.04 

0.68 

0.08 

1.0 

0.104 

0.094 

0.11 

425 

0.14 

0.04 

0.68 

0.08 

1.5 

0.135 

0.105 

0.12 

417 

0.139 

0.04 

0.68 

0.079 

2.45 

0.148 

0.12 

0.118 

362 

0.14 

0.04 

0.68 

0.08 

Electrolyte  thickness  varied  (anode  support  thickness  = 
Electrolyte  thickness  (|xm) 

1 .0  mm,  anode  support  porosity  = 

48%,  cathode  interlayer  thickness 

=  20  |j.m) 

4 

0.103 

0.09 

0.13 

500 

0.14 

0.04 

0.68 

0.08 

8 

0.104 

0.094 

0.11 

425 

0.14 

0.04 

0.68 

0.08 

15 

0.114 

0.088 

0.095 

396 

0.14 

0.04 

0.68 

0.08 

20 

0.142 

0.086 

0.095 

404 

0.14 

0.042 

0.68 

0.082 

Anode  support  porosity  varied  (Electrolyte  thickness  =  8  |xm,  anode  support  thickness  =1.0  mm,  cathode  interlayer  thickness  =  20  p,m) 


Anode  support  porosity 


32% 

0.149 

0.096 

0.092 

352 

0.14 

0.04 

0.22 

0.085 

48% 

0.104 

0.094 

0.11 

425 

0.14 

0.04 

0.68 

0.08 

57% 

0.094 

0.03 

0.079 

684 

0.14 

0.04 

0.75 

0.08 

76% 

0.08 

0.022 

0.09 

783 

0.14 

0.04 

0.82 

0.08 

Cathode  interlayer  thickness  varied  (electrolyte  thickness  =  8  |xm,  anode  support  thickness  =  1.0  mm,  anode  support  porosity  =  48%) 
Cathode  interlayer  thickness  (|xm) 


6 

0.096 

0.176 

0.155 

321 

0.135 

0.04 

0.61 

0.08 

20 

0.104 

0.094 

0.11 

425 

0.14 

0.04 

0.68 

0.08 

56 

0.118 

0.095 

0.155 

543 

0.135 

0.04 

0.61 

0.08 

105 

0.134 

0.105 

0.155 

598 

0.138 

0.04 

0.63 

0.08 

of  (~45%  porosity)  and  £>02-n2  (~26%  porosity) 

from  curve-fitting  as  ~0. 14  and  ~0.04cm2  s'1,  respectively, 
are  in  excellent  agreement  with  the  actual  measurements. 
Also  note  that  .D^-h20  >  ,  consistent  with  expecta¬ 
tions  since  >  7?o2-N2  •  by  virtue  of  the  much  lower 

molecular  weight  of  Fb  as  compared  to  the  other  species.  At 
800°C,  the  Z>h2-h2o  is  ~7.7cm2s“1,  estimated  using  the 
Chapman-Enskogg  model  [14].  The  anode  support  porosity 
is  ~48%.  Based  on  the  estimated  Onf-HiO  by  curve  fitting, 
the  measured  porosity,  and  estimated  Z>h2-h2o>  a  totuosity 
factor  of  ~5.4  is  determined,  which  is  quite  reasonable. 

Also  listed  in  Table  5  are  the  corresponding  a  and  b, 
which  describe  the  net  activation  polarization,  wherein  the 
activation  polarizations  at  both  the  anode  and  the  cathode  are 
lumped  into  one.  As  the  activation  polarizations  from  both 
electrodes  are  lumped  into  one  single  polarization  equation, 
the  parameters  a  and  b  are  phenomenological  fitting  parame¬ 
ters,  and  the  exchange  current  density  obtained  from  these  is 
also  a  phenomenological  fitting  parameter.  These  parameters 
using  the  Butler- Volmer  type  equation  can  be  given  by 


a 


RT 
4  aF 


In  i0 


and 

RT 

b  &  - 

4  a  F 


(21) 


(22) 


where  a  is  the  transfer  coefficient.  Thus,  from  Eqs.  (21)  and 
(22),  the  exchange  current  density  is  given  in  terms  of  the 


fitting  parameters  a  and  b  by 

i°  =  Btp  (~l) 


(23) 


It  is  important  to  note  that  for  the  fitting  procedure  and 
units  used,  there  will  be  corresponding  units  for  i0,  here 
A  cm'2.  That  is,  the  parameters  a  and  b  depend  upon  the 
choice  of  units,  and  thus  so  do  the  units  of  i0.  Thus,  for 
example,  for  an  anode  support  thickness  of  0.5  mm  with 
a  =  0.098  and  ft  =  0.11,  the  i0  =  (-0.098/0. 1 1)  =  0.41  A  cm'2 
or  410mA  cm'2.  Table  5  shows  that  the  estimated  ex¬ 
change  current  densities  are  essentially  the  same,  about 
415  mAcm'2  virtually  for  all  anode  support  thicknesses,  ex¬ 
cept  for  anode  support  thickness  of  2.45  mm,  for  which  it  is 
somewhat  lower  (~362mAcm“2).  The  significance  of  the 
lower  ia  for  thick  anode  supports  will  be  discussed  later.  Ex¬ 
amination  of  Fig.  18  shows  that  fitted  curve  diverges  below 
about  0.3  A  cm'2,  consistent  with  the  expectation  that  for 
current  densities  lower  than  about  i0,  the  Tafel  equation  does 
not  satisfactorily  describe  activation  polarization. 

Table  5  also  gives  the  fitting  parameters  for  set  of  cells 
tested  at  800  °C  wherein  the  only  parameter  varied  was  the 
electrolyte  thickness,  all  other  parameters  maintained  at  the 
standard  values.  The  values  of  R,  obtained  from  the  best  fits 
vary  between  0. 103  and  0. 142  Q  cm2  for  an  electrolyte  thick¬ 
ness  ranging  between  4  and  20  |xm.  These  compare  well  with 
the  ohmic  ASR  determined  by  current  interruption  given  in 
Table  1,  which  ranges  between  0.1  and  0.14  Q  cm2.  The  en¬ 
tire  set  of  data  could  be  fitted  well  with  the  same  values  of 

neff(l)  r-effra)  r.eff(l)  ,  r.eff(2)  ,  .  ,  . 

£>02-N2’  £>o2-N2’  DH2-H20’  and  DH2-H20’  Whlch  1S  consis- 
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tent  with  expectations  since  the  microstructures  and  porosi¬ 
ties  of  the  anode  support,  anode  interlayer,  cathode  current 
collector  and  cathode  interlayer  were  the  same  in  all  of  these 
cells.  The  fitted  parameters  a  and  b  give  the  magnitude  of  the 
exchange  current  density,  i0,  which  ranges  between  ~400  and 
~500mAcm-2,  the  same  as  before. 

Table  5  further  gives  the  results  of  fitting  data  on  cells 
wherein  the  only  parameter  varied  was  the  anode  support 
porosity,  all  other  parameters  fixed  at  the  standard  values. 
The  R{  values  corresponding  to  the  best  fits  vary  between 
0.149  cm2  for  anode  support  porosity  of  32%  to  0.08  !Tcm2 
for  anode  support  porosity  of  76%.  As  shown  in  Table  1,  the 
ohmic  ASR  obtained  using  current  interruption  for  the  same 
porosity  range  varies  between  0.156  and  0.074  £2  cm2,  indi¬ 
cating  good  agreement.  The  transport  parameters  correspond¬ 
ing  to  the  cathode  current  collector,  ,  cathode  inter¬ 

layer,  Z?q^n2,  and  anode  interlayer,  77h*-h20'  aie  essen¬ 
tially  the  same  as  before,  consistent  with  expectations  since 
the  respective  microstructures  of  these  regions  were  the  same 
in  all  cells  tested.  Note,  however,  that  the  which 

is  the  effective  diffusivity  through  the  anode  support,  varies 
between  0.22  cm2  s_  1  for  a  cell  with  anode  support  porosity 
of  32%  to  0.82  cm2  s-1  for  anode  support  porosity  of  76%. 
The  profound  effect  of  anode  support  porosity  on  H2-H2O 
gaseous  transport  and  the  corresponding  anodic  concentra¬ 
tion  polarization  is  evident.  The  parameters,  which  describe 
activation  polarization,  namely  a  and  b,  are  also  listed  in 
Table  5.  It  is  interesting  to  note  that  the  estimated  exchange 
current  density,  i0,  varies  between  ~352mAcm-2  for  an¬ 
ode  support  porosity  of  32%,  to  ~783mAcm“2  for  anode 
support  porosity  of  76%.  This  result  leads  to  the  following 
observations/conclusions.  (1)  The  fact  that  changes  made  on 
the  anode  side  affect  the  overall  exchange  current  density,  i0, 
suggests  that  there  must  be  significant  contribution  from  both 
the  cathode  and  the  anode  to  the  net  activation  polarization. 
(2)  The  anode  interlayer,  where  most  of  the  anodic  electro¬ 
chemical  reaction  occurs,  was  the  same  in  all  cells.  Yet  the 
anodic  activation  polarization  appears  to  be  affected.  This 
apparent  conflict  can  be  resolved  when  the  factors  which  in¬ 
fluence  activation  polarization  are  examined  more  closely. 
It  is  to  be  noted  that  activation  polarization  is  a  function 
of  the  concentration  (partial  pressure)  of  the  electro-active 
species  in  the  region  of  the  electrode  where  electrochemical 
reactions  occur  (i.e.,  in  the  interlayer);  the  higher  the  con¬ 
centration  (partial  pressure)  of  the  electro-active  species,  the 
higher  should  be  the  exchange  current  density,  iQ,  and  the 
lower  should  be  the  activation  polarization.  For  a  given  im¬ 
posed  current  density,  the  lower  is  the  anode  support  porosity, 
the  lower  will  be  pH2(i)(f)  and  p\i2^(i),  the  lower  will  be  the 
exchange  current  density,  iQ,  and  the  higher  will  be  the  acti¬ 
vation  polarization.  That  is,  the  activation  polarization  is  not 
independent  of  concentration  polarization.  For  this  reason,  as 
the  anode  support  porosity  is  decreased,  there  also  is  an  effect 
on  the  estimated  i0.  The  observation  that  for  a  cell  with  anode 
support  thickness  of  2.45  mm  the  exchange  current  density 


was  lower  than  for  thinner  anode  supports,  is  also  consistent 
with  this  reasoning.  The  implication  is  that  Eq.  (18),  which 
has  been  used  to  describe  the  overall  V(i)  versus  i  polarization 
behavior  is  rather  crude,  since  it  ignores  the  dependence  of  a 
and  b  on  the  partial  pressures  of  the  electro-active  species;  for 
the  anode  side  pn2(i)(i)  and  P'w-AiFb  and  for  the  cathode  side 
po2(i)(f)  and  p'q2(i\(j)-  This  also  implies  that  more  complex 
models  which  incorporate  the  complete  Butler- Volmer  equa¬ 
tion,  such  as  in  the  work  of  Zhu  and  Kee,  may  be  required  to 
demonstrate  better  quantitative  agreement  [10]. 

Table  5  gives  the  results  of  curve  fitting  for  cell  tests 
wherein  the  cathode  interlayer  thickness  was  varied  between 
~6  and  ~105  |xm,  all  other  parameters  fixed  at  the  standard 
values.  Once  again  the  values  of  -n2  -  °o!' -n2  >  T^-HzO’ 

and  T>h^h->0  corresponding  to  the  best  fits  are  essentially 
the  same  as  for  all  other  cells  with  the  same  microstructures 
and  porosities  (that  is  the  standard  cells,  those  with  different 
electrolyte  thicknesses,  and  cells  with  different  anode  sup¬ 
port  thicknesses).  The  ohmic  ASR  values  from  the  data  fit  are 
also  consistent  with  those  determined  by  current  interruption. 
From  the  fitted  parameters  a  and  b,  the  estimated  exchange 
current  density,  i0,  for  cathode  interlayer  thickness  of  ~6  |xm 
is  about  321  mAcm-2,  which  is  the  lowest.  For  cathode  in¬ 
terlayer  thickness  of  20  pm,  the  iQ  is  ~425  mA  cm-2,  while  it 
increases  to  over  500  mA  cm-2  for  cathode  interlayer  thick¬ 
nesses  of  56  and  105  pm.  Theoretical  analysis  by  Tanner  et 
al.  [15]  for  a  composite  electrode  shows  that  effective  charge 
transfer  resistance  is  given  by 


Rct(eff) 


PidR° t 
1  -  Vv 


(24) 


where  R°t  is  the  intrinsic  charge  transfer  resistance  (which  is 
representative  of  the  electrocatalyst  spread  evenly  on  a  planar 
electrolyte  surface),  p\  is  the  ionic  resistivity  of  the  composite 
electrode  and  d  is  the  grain  size.  The  measured  effective 
exchange  current  density,  alternatively,  can  be  given  by 


RT  I  1  -  Vy 
4  F]j  pi  R°td 


(25) 


Both  of  these  two  Eqs.  (24)  and  (25),  are  asymptotic  limits 
of  the  general  equations  given  by  Tanner  et  al.  [15],  which 
shows  that  as  the  electrode  interlayer  thickness,  h,  increases, 
the  exchange  current  density  increases  (or  effective  charge 
transfer  resistance  decreases),  and  for  h  >  I0d-I5d  reaches 
an  asymptotic  value  given  by  Eq.  (25)  (or  Eq.  (24)).  In  the 
present  experiments,  the  typical  grain  size  in  the  cathode  in¬ 
terlayer  was  ~2  pm.  Hence,  it  is  expected  that  for  an  inter¬ 
layer  thickness  above  about  20-30  pm,  the  asymptotic  value 
would  be  achieved.  The  observation  that  the  exchange  current 
density  increases  with  increasing  cathode  interlayer  thick¬ 
ness  is  consistent  with  this  model.  Note  that  the  i0  seems  to 
reach  an  asymptotic  value  for  cathode  interlayer  thickness  of 
~50  pm.  An  interesting  observation  is  that  even  though  the 
cell  with  20  pm  cathode  interlayer  has  a  lower  i„  than  the  cell 
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with  ~56  |JLm  cathode  interlayer,  it  exhibits  higher  power  den¬ 
sity.  This  is  consistent  with  the  fact  that  the  cell  with  56  |xm 
thick  cathode  interlayer  exhibits  higher  ohmic  ASR,  which 
seems  to  arise  from  the  cathode  interlayer.  In  addition,  the 
larger  the  cathode  interlayer  thickness,  the  higher  is  the  ca¬ 
thodic  concentration  polarization.  The  model  by  Tanner  et 
al.  does  not  account  for  ohmic  losses  in  the  cathode  inter¬ 
layer,  which  can  occur  due  to  inhomogeneous  distribution  of 
particles  present  in  real  composite  electrodes  [15]. 

4.3.  Temperature  dependence  of  polarization 

The  ohmic  polarization  is  temperature  dependent  mainly 
by  virtue  of  the  thermally  activated  dependence  of  the  YSZ 
ionic  resistivity;  the  lower  the  temperature,  the  higher  is  the 
resistivity.  Note  that  the  ASR  for  the  standard  cells  at  800  °C 
is  0. 104  £2  cm2  while  that  at  700°  C  is  0. 19  cm2.  The  acti¬ 
vation  polarization  is  also  thermally  activated;  which  is  re¬ 
flected  in  the  thermally  activated  dependence  of  the  exchange 
current  density  i0  through  two  temperature-dependent  pa¬ 
rameters,  namely,  p\  and  R°v  The  concentration  polarization, 
however,  is  weakly  dependent  on  temperature.  There  are  two 
sources  of  temperature  dependence  of  concentration  polariza¬ 
tion.  One  of  them  is  the  temperature  dependence  of  the  partial 
pressure  of  hydrogen  at  the  anode  interlayer/electrolyte  inter¬ 
face  (Eq.  (8)),  and  partial  pressure  of  oxygen  at  the 

cathode  interlayer/electrolyte  interface  (Eq.  (13)),  p'^^ii), 
which  exhibit  temperature  dependence  directly  as  well  as 
through  the  temperature  dependence  of  effective  diffusivi- 
ties.  The  binary  diffusivities,  through  the  Chapman-Enskogg 
equation  are  proportional  to  T312 .  The  other  dependence  is  di¬ 
rect  linear  relationship  with  temperature,  as  reflected  in  Eq. 
(7)  for  >7a,con  and  Eq.  (18)  for  rjcxon.  The  temperature  depen¬ 
dence  of  concentration  polarization  thus  is  generally  weak. 
Calculations  of  ^a,con  and  ()CiCOn  made  using  Eqs.  (7)  and  (18), 
respectively,  show  that  these  are  weakly  dependent  on  tem¬ 
perature  over  the  temperature  range  between  600  and  800  °C. 
Fig.  19,  for  example,  shows  calculated  );a,con  as  a  function  of 
current  density  over  a  range  of  temperatures  between  600 
and  800  °C.  Note  that  the  ?;acon  slightly  increases  with  tem¬ 
perature.  Thus,  the  principal  temperature  dependence  of  cell 
performance  is  due  to  the  temperature  dependence  of  ohmic 
and  activation  polarizations. 

4.4.  Estimation  of  anode  and  cathode-limiting  current 
densities  and  concentration  polarization  as  a  function  of 
various  cell  parameters 

From  the  estimated  Dh!-h20’  and 

^H2-H20  from  curve-fitting  to  V(i)  versus  i  data,  the  zas  and  ics 
were  estimated  using  Eqs.  (10)  and  ( 14),  respectively,  for  cells 
with  different  anode  support  thicknesses  and  anode  support 
porosities.  The  estimated  zas  and  z'cs  are  listed  in  Table  6.  Us¬ 
ing  the  anc*  ^[11-iioO'  anode  concentration  polar¬ 

ization  was  calculated  for  a  number  of  cases.  Also  calculated 


Fig.  19.  Anode  concentration  polarization  for  the  cell  with  optimized  pa¬ 
rameters  (0.5  mm  anode  support  thickness,  57%  anode  support  porosity)  as 
a  function  of  current  density  over  a  temperature  range  between  600  and 
800  °C. 

was  the  partial  pressure  of  hydrogen  at  the  electrolyte/anode 
interlayer  interface,  Plots  of  the  calculated  /?jh(()(z) 

and  /;a, con  versus  i  are  given  in  Fig.  20  for  anode  support 
thickness  ranging  between  0.5  and  2.45  mm.  The  profound 
effect  of  anode  support  thickness  on  anodic  concentration 
polarization  is  clearly  seen  in  the  figure.  Thus,  making  the 
anode  support  as  thin  as  possible  is  important,  although  from 
the  standpoint  of  mechanical  ruggedness,  probably  the  lower 
limit  is  about  0.5  mm.  It  is  also  to  be  noted  that  the  decrease  in 
P\i,t  j){i)  at  higher  anode  support  thicknesses  leads  to  higher 
activation  polarization  since  a  and  b  are  functions  of  partial 
pressures  of  hydrogen  (at  the  anode  interlayer/electrolyte  in¬ 
terface).  Fig.  21  shows  similar  plots  of  p'^2fi)(i)  and  z;aiCon 
versus  i  where  the  parameter  varied  was  the  anode  support 
porosity.  Note  the  large  effect  of  anode  support  porosity,  es¬ 
pecially  at  low  values  of  porosity.  In  Figs.  20  and  21,  the 
vertical  scale  on  the  left  axis  is  in  mV.  Note  that  concentra¬ 
tion  polarization  as  high  as  ~  150  mV  (or  even  greater)  is  esti¬ 
mated  at  high  current  densities.  Fig.  22  shows  similar  plots  of 
^2(o(0  and  7c, con  versus  i  for  the  case  where  the  only  param- 


Table  6 

Limiting  current  densities  for  various  anode  thicknesses  and  anode  porosities 


ias  (A  cm  2) 

zcs  (A  cm  2) 

Anode  thickness  (mm) 

0.5 

21.81 

12.8 

1.0 

12.49 

1.5 

8.75 

2.45 

5.58 

Anode  porosity  (%) 

32 

4.49 

48 

12.49 

57 

13.57 

76 

14.63 

94 
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Fig.  20.  Anode  concentration  polarization,  r]^conc,  and  partial  pressure  of 
hydrogen  at  the  anode  interlayer/electrolyte  interface,  Ph2(,-)(0»  as  a  function 
of  current  density  for  cells  with  different  anode  support  thicknesses. 

eter  varied  was  the  cathode  interlayer  thickness.  In  general, 
the  cathodic  concentration  polarization  is  small  as  long  as 
the  cathode  interlayer  thickness  is  small.  However,  for  large 
cathode  interlayer  thicknesses  (e.g.  —50-100  |xm),  the  r/c  con 
was  estimated  to  be  as  high  as  —20  mV.  Especially  significant 
is  the  fact  that  the  corresponding  p'(h(jXi)  can  be  substantially 
lower  than  ,  which  thus  also  increases  activation  polar¬ 
ization,  since  a  and  b  depend  upon  on  Pq2 (,)(*)■  Thus,  even  in 
an  anode-supported  cell,  it  is  important  to  ensure  that  oxidant 
transport  through  the  cathode  is  not  impeded. 

4.5.  The  optimized  cell 

Using  the  maximum  power  density  (MPD)  as  a  figure 
of  merit,  note  that  the  optimized  cell  exhibited  an  MPD 


Current  Density  (A/cm2) 


Fig.  21.  Anode  concentration  polarization  t]a,conc.  and  partial  pressure  of 
hydrogen  at  the  anode  interlayer/electrolyte  interface,  pjj2(;)0). as  a  function 
of  current  density  for  cells  with  different  anode  support  porosities. 
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Fig.  22.  Cathode  concentration  polarization,  r/c  conc,  and  the  partial  pressure 
of  oxygen  at  cathode  interlayer/electrolyte  interface,  Po2(i)0)’  as  a  function 
of  cell  current  density  for  cells  with  different  cathode  interlayer  thicknesses. 

of  ~1.8  W  cm-2  at  800  °C,  ~0.8  W  cm”2  at  700  °C  and 
— 0.4  W  cm-2  at  600  °C.  All  three  polarization  losses  were 
proportionately  lower  in  the  optimized  cell  as  compared  to 
the  standard  cell  which  exhibited  MPD  of  —  1.2Wcm-2 
at  800  °C,  —0.58  W  cm-2  at  700  °C  and  -0.2  W  cm”2  at 
600  °C.  In  the  present  work,  only  few  parameters  were  varied, 
namely  the  various  thicknesses,  and  anode  support  porosity. 
It  is  clear  that  further  performance  gains  are  to  be  expected 
when  the  dependence  of  various  polarization  losses  on  other 
parameters  is  examined,  such  as  the  relative  proportions  of 
the  phases  in  the  composite  electrodes,  and  microstructures 
of  the  electrodes.  The  observation  that  even  within  a  given 
materials  set  the  cell  performance  can  be  varied  over  a  wide 
range  exemplifies  the  role  of  various  parameters  on  polar¬ 
izations,  and  underscores  the  importance  of  investigating  the 
role  of  microstructural  and  geometric  parameters.  It  also  sug¬ 
gests  that  when  investigating  the  effect  of  various  materials 
on  cell  performance,  a  comparison  among  different  materi¬ 
als  by  itself  may  be  of  little  value,  unless  a  careful  study  of 
the  microstructural  and  geometric  parameters  is  also  made 
concurrently. 


5.  Summary 

In  the  present  work,  a  series  of  anode-supported  cells  com¬ 
prising  a  mixture  of  Ni  +  YSZ  anode  support,  Ni  +  YSZ  an¬ 
ode  interlayer,  YSZ  electrolyte,  LSM  +  YSZ  cathode  inter¬ 
layer,  and  LSM  cathode  current  collector,  were  fabricated 
with  one  parameter  varied  at  a  time.  The  cell  parameters  in¬ 
vestigated  were  the  electrolyte  thickness,  the  anode  support 
thickness,  the  cathode  interlayer  thickness,  and  anode  sup¬ 
port  porosity.  Electrochemical  performance  was  measured 
with  humidified  hydrogen  as  fuel  and  air  as  oxidant  over 
a  temperature  range  from  600  to  800  °C,  with  most  of  the 
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testing  done  at  800  °C.  The  ohmic  ASR  was  measured  by 
current  interruption.  It  was  observed  that  the  various  polar¬ 
izations  and  ohmic  ASR  are  functions  of  the  cell  parameters 
investigated.  The  electrochemical  performance  of  cells  was 
fitted  to  a  simple  analytical  model.  The  emphasis  was  on  de¬ 
scribing  voltage  versus  current  density  performance  curves 
with  the  minimum  possible  number  of  parameters,  which  in 
principle  can  be  independently  measured.  The  results  showed 
that  a  significant  contribution  to  the  ohmic  ASR  is  attributed 
to  components  other  than  the  electrolyte  at  temperatures  as 
low  as  700  °C,  as  long  as  the  electrolyte  thickness  is  small 
(~  1 0  pan).  This  shows  that  YSZ  as  electrolyte  is  satisfac¬ 
tory  in  a  thin  film  form  down  to  at  least  700  °C,  and  perhaps 
down  to  600  °C.  Through  the  optimization  study,  a  set  of  pa¬ 
rameters  was  identified,  an  optimized  cell  was  fabricated  and 
its  electrochemical  performance  was  evaluated.  The  maxi¬ 
mum  power  density  of  the  optimized  cell  was  ~1.8W  cm-2 
at  800  °C  and  ~0.4  W  cm-2  at  600  °C.  The  results  show  that 
a  significant  improvement  in  cell  performance  can  be  real¬ 
ized  by  addressing  microstructural  and  geometric  aspects  of 
cell  parameters,  even  when  using  the  standard  set  of  mate¬ 
rials  (LSM,  YSZ,  and  Ni).  Further  performance  gains  are 
to  be  expected  with  optimization  of  other  cell  parameters, 
which  include  the  anode  interlayer  thickness,  cathode  cur¬ 
rent  collectors  and  also  the  relative  proportions  of  phases 
in  the  composite  electrodes.  It  is  clear  that  even  greater  im¬ 
provements  are  possible  using  materials  with  intrinsically  su¬ 
perior  properties.  The  present  study  underscores  the  impor¬ 
tance  of  microstructural  and  geometric  parameters,  and  the 
need  to  carry  out  a  quantitative  optimization  study  for  a  given 
materials  set. 
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